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A strong light pulse propagating in a nonlinear Kerr medium produces a change in the refractive 
index, which makes light travel at different speeds inside and outside the pulse. By tuning the pulse 
velocity, an analogue black hole horizon can be obtained in a suitable frequency window. In this 
paper, we develop a quantum theory of light propagation for this system, including the frequency 
dispersion of the refractive index of the medium by coupling the electromagnetic field to matter 
polarization fields. In a configuration with a single black horizon, the spectrum of spontaneously 
emitted particles presents a few similarities with Hawking radiation. In horizonless system sponta- 
neous vacuum emission is still possible due to the dispersive nature of the medium, yet with different 
spectral properties. 
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I. INTRODUCTION 

Hawking radiation [1, 2] is the quantum production 
of particles from vacuum fluctuations due to the pres- 
ence of a black hole horizon in a curved stationary ge- 
ometry. In the pioneering work [3], it was shown that 
this prediction does not reside on the peculiar dynamical 
features of a general-relativistic spacetime, but only on 
the kinematic properties of a quantum field living in a 
curved spacetime. As a result, analogous quantum vac- 
uum emission phenomena have been anticipated to occur 
in several physical systems, ranging from flowing fluids 
and superfluids, to ion rings, to nonlinear optical sys- 
tems [4]. In particular, the advanced techniques of pulse 
manipulation and light detection that have been devel- 
oped in the recent years put nonlinear optical systems 
among the most promising candidates for the realization 
of analogue models of gravitational systems. 

In addition to the many proposals that appeared in 
the last few years [5-8], the first claim of observation 
of an analogue Hawking radiation in a laboratory was 
reported in [9, 10] using the refractive index change in- 
duced by a strong laser pulse propagating across a non- 
linear dielectric medium. The velocity of the pulse can 
be tuned either by changing the laser wavelength or by 
using an axicon with different angles to produce a Bessel 
beam [11-13]). In this way the speed of optical frequency 
photons inside the pulse can be made smaller than the 
pulse velocity, while the one outside the pulse remains 
larger. Consequently, the boundaries of the pulse ap- 
pear as analogue horizons, as seen from the frame co- 
moving with the pulse. Unfortunately, the observation 
of Hawking emission in this system is still considered as 
controversial by some authors, who recently raised a few 
issues [14, 15]. Inspite of attemps of alternative theoreti- 
cal interpretation [16, 17], no satisfactory explanation of 
the experimental observations has been found yet. 

Differently from previous studies that deeply reside 
on the analogy with gravitational physics, the present 



work develops a microscopical quantum optical model of 
light propagation in a nonlinear dielectric modulated by 
the passage of a high-intensity laser pulse. The struc- 
ture of modes propagating in a one-dimensional system 
is described by using the Lagrangian of the electromag- 
netic field coupled with three polarization fields. In this 
way, we are able to canonically quantize a theory which 
exactly reproduces the complete Sellmeier dispersion of 
transparent materials such as the fused silica used in the 
experiment [9] . This allows us to describe the system in a 
far more realistic way than existing works based on sim- 
plified subluminal dispersions [18, 19] and obtain quan- 
titative predictions for the spectrum of spontaneously 
emitted photons under the simplifying assumption that 
the strong pulse propagates through the medium in a 
steady and rigid way and that the jump in the spatial 
profile of the refractive index is very sharp. 

In particular, we point out how the highly nontriv- 
ial dispersion relation makes the analogy with standard 
quantum field theory in curved spacetime much weaker 
than in other analogue systems considered in the lit- 
erature, such as Bose-Einstein condensates [4]: for in- 
stance, in the present nonlinear optical context, an ana- 
logue horizon can be defined only in a finite frequency 
range [20] that does not extend to the long-wavelength 
limit where the analogue geometry is generally defined. 
Nonetheless, the properties of the quantum vacuum emis- 
sion in this frequency range share several features of stan- 
dard Hawking radiation. 

Furthermore, in contrast to non-dispersive media, 
where a necessary and sufficient condition to trigger emis- 
sion processes a la Hawking is the presence of analogue 
horizons, quantum vacuum radiation occurs in dispersive 
media even in the absence of any horizon as soon as the 
dispersion allows for modes with a negative norm (with 
respect to the natural scalar product induced by the 
symplectic structure of the Hamiltonian) and a positive 
frequency as seen from the pulse comoving frame [20]. 
With respect to the above mentioned Hawking-like emis- 
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sion channel, this additional emission is generally much 
weaker. 

The structure of the article is the following. In Sec. II, 
we start from the Lagrangian to derive and solve the 
equation of motion for the electromagnetic field coupled 
to matter polarization fields. In Sec. Ill, the properties 
of the eigenmodes are investigated in the most significant 
configuration with a single analogue black hole horizon, 
and the two bases of in- and out-going scattering modes 
are built. In Sec. IV, the spectrum of emitted particles 
is computed in the case of a large refractive index jump 
between the two sides of the analogue black hole horizon. 
The more realistic case of a small refractive index jump 
is studied in Sec. V along with the horizonless configu- 
ration recently realized the experiment [9]. Conclusions 
are finally drawn in Sec. VI. 



II. THE GENERAL FRAMEWORK 



In this section we develop the quantum theory used to 
describe light propagation in a dielectric medium with 
a refractive index modulation that moves in space at a 
uniform velocity v. in the experiments, the refractive 
index change is generated by the strong pulse via the 
Kerr nonlinearity of the medium. At a simplest level of 
approximation, one can assume that the refractive index 
change follows locally in space and instantaneously in 
time the intensity profile of the pump pulse [21]. 

The optical properties of the medium are described in 
the spirit of the Hopfield model [22] in terms of an elec- 
tromagnetic field interacting with a polarization field. To 
closely reproduce the Sellmeier dispersion relation of typ- 
ical transparent dielectrics such as fused silica [9], the 
polarization of the medium has to show three poles with 
different strengths and frequencies. The effect of the 
passage of the pulse in the medium is then modeled as 
a spatio-temporal variation of these quantities, consid- 
ered as external parameters. More sophisticated models 
where the internal level structure of the emitters is ex- 
plicitly taken into account were developed in [23-25] but 
the complexity of such a description goes far beyond the 
scope of the present work. 

For the sake of simplicity, the strong pulse is assumed 
to propagate across the medium at a constant speed 
and to maintain a constant shape in space. We also 
restrict our description to a one-dimensional geometry 
where light modes propagate parallel to the direction of 
the pulse velocity. The refractive index jump is assumed 
to be localized at the sharp interface between two asymp- 
totic homogeneous regions. More general configurations 
with smooth refractive index profiles will be the subject 
of future work. 




A. The field equation 

In the laboratory reference frame, the Lagrangian den- 
sity in one dimension of the electromagnetic field coupled 
to N polarization fields Pi is 



(l/ 

where A and Pi oscillate in a direction orthogonal to 
their propagation direction. In this simple model the 
polarization is described by fields of harmonic oscillators, 
with strength /3 _1 and inertia (Piilj) -1 . For the sake 
of simplicity, in this paper we shall restrict to the case 
N = 3, which suitably describes the dispersion relation 
in fused silica, the material used in the experiment of 
Ref. [9]. 

A propagating strong light pulse causes a local per- 
turbations of the parameters /3j and f2j. Consequently, 
the system is stationary when observed in the reference 
frame comoving with the pulse at velocity v. Thus, it 
is convenient to transform the laboratory coordinates X 
and T to the comoving coordinates x and t, by applying 
a Lorentz boost A with velocity v 



t = 1 [T-vX/c% 
x = j[X — vT], 



(2) 



where 7 = 1/ \Jl — v 2 /c 2 . Coherently, differential opera- 
tors transform as 



d T =j[d t - vd x ], 
d x = 1 [d x -vd t /c 2 }, 



(3) 



Furthermore, we shall treat A and Pi as scalar fields, i.e., 
we do not transform them under the boost. Although this 
might appear not correct, it is completely legitimate, as 
shown in Appendix A. 

The transformed Lagrangian density is 



L = 



A 2 



A' 



3 r 



8ttc 2 8tt 
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P? 



Pi - vP' 



2/3, c 



-A 



, (4) 



where dot and prime denote derivation with respect to t 
and x, respectively. 

As usual, the conjugate momenta are obtained by vary- 
ing the Lagrangian 



L = / dx L 



(5) 



with respect to the time derivatives of A and Pf. 



n A = 



A 

47TC 2 ' 



Hp, = 



P l -vP' i )+^A. (6) 
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We can impose canonical commutation relations on A 
and Pi (see Appendix A 



[A(x),U A (x')] = ih6{x-x / ), 
[P i {x),Hp j (x')} = ih6 ij 5(x-x') 

and all the other commutators vanish. 
The Hamiltonian density 



n = 



1 



3 

A u A + u A A + n p s + n P A 



becomes 



(7) 



-C (8) 



where I± is the 4x4 identity matrix. With this notation 
the Hamilton equations can be written in a compact form 
as 



V = 7?(WH). 

Moreover, a scalar product 



(17) 



(Vi,v 2 : 



dxV;{x,t)r]V 2 (x,t) 



(18) 

can be defined on the set of the solutions of Eq. (17), 
generalized to complex values. In virtue of Eq. (17) , using 
rf = —Is, t^t] — 7g, and the fact that the Hamiltonian 
density "H is quadratic in the fields and their momenta, 
this scalar product is conserved by time evolution: 



H = 2ttc 2 nl 
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n, 
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-(p/n Pi + n P( ^) 



(9) 



and the Hamilton equations are derived by the commu- 
tators of the fields and their conjugate momenta with the 
Hamiltonian 



H = / dxU. 



We obtain 

A = 4ttc 2 n^, 



n. 



n, 



A" 
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^4 + vP', 
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(12) 

(13) 

(14) 



Note that one might have obtained the first two equations 
directly from the definition of conjugate momenta (6). 

It is now convenient to define the eight-dimensional 
vector 

V = (A P x P 2 P 3 U A U Pl Up 2 U Pa ) T (15) 
and the matrix 



V 



h 

-h o 



(16) 



d t (V 1 ,V 2 
dx 



d t (V;)r ] V 2 + V?r ] d t (V 2 ) 



(y v n)\ v=Vi) v 2 -v? (y v n)\ v=V2) 



= o, 

(19) 



Being the system stationary in the reference system 
comoving with the pulse, the real field V can be expanded 
on a basis of frequency eigenmodes V^: 



V = 



Jo 



where 



(20) 



(21) 



the label a denotes various modes with the same eigen- 
frequency uj, and V" are properly normalized (see Ap- 
pendix B) with respect to the norm induced by the scalar 
product defined in Eq. (18). 



B. Homogeneous systems 

In the asymptotic regions, far from the perturbation, 
the system is homogeneous and the parameters f2j, /3j, 
and v are constant both in time and space. In this situ- 
ation one can chose V™ as momentum eigenmodes: 

V°{x,t) = e-^t+^xya^ ^ 

where V" is a vector of constant C-numbers, satisfying 
-iuV£ =T)lC(k a )V2, (23) 

and 
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(24) 



Equation (23) can be written as 

{ir,K(k a )-cjI 8 )V« = 0. (25) 
which has nonvanishing solutions if and only if 

det(ir]JC(k a ) - = 0. (26) 

The computation of the above determinant yields the dis- 
persion relation 



,2;„2 



c z k 



A Airfc^iuj + vk) 2 

^ i--f 2 (uj + vk) 2 /n 2 ' 



(27) 



that has in general 8 solutions fc Q for each value of the 
frequency u. 

Applying the boost of Eq. (2) to to and fc, the corre- 
sponding frequency and wavenumber K, measured in 
the laboratory reference frame, are 



j[oj + vk], 
j[k + vlu/c 2 ]. 



(28) 



Equation (27) becomes 

c 2 k 2 = n 2 



4ttA 



n 2 /n 2 



(29) 



which is the well-known Sellmeier dispersion relation [26, 
27]. 

The solutions of the system (23) (eigenmodes of ir/K.) 
are 



/ c \ 

ipx 7(w + vk a )[l - 7 2 (w + vkf/n 2 ]- 1 

ifa 7 (oj + »fc a )[l - 7 2 (w + wfc) 2 /^]-! 
7 (w + wfc a )[l - 7 2 (^ + ^fc) 2 /"!]" 1 
— ius /Aire 
7[l- 7 2 (w + wfc) 2 /fi 2 ]- 1 
7 [l-7 2 (w + «fc) 2 /0 2 ]- 1 
7[l- 7 2 (w + wfc) 2 /0 2 ]- 1 



K Q = C°. 



V 

where the constant C" is computed in Appendix B 

47r/3i7 2 (o; - wfc„) 



(30) 



C fc a 



212 



(31) 



using the normalization condition 

I teC, 2 ) I = ^2 ~ 



(32) 



The modulus is needed since some modes have negative 
norm, because the scalar product of Eq. (18) is not pos- 
itive definite. In particular, it is possible to show (see 
Appendix B) that the norm of a mode is positive when 
its laboratory frequency O is positive, whereas it is nega- 
tive when £1 < 0, independently of the value of K. Since, 
for a given positive value of the comoving frequency u>, 
the dispersion relation (27) admits solutions for k, such 
that n = "f(uj + vk) is negative, there exist modes with 
positive frequency w, but negative norm. This implies 
that, in the expansion of the field V, the Fock operators 
associated with those positive-o; modes are not destruc- 
tion but instead creation operators. Naming P the set of 
positive norm modes V", labeled by a, and N the set of 
negative norm modes V", labeled them by a, V becomes 



V 



du e 




- e +ik ' ■■' 



+ h.c, (33) 



where h.c. stands for hermitian conjugate. Note that the 
positive frequency part of the field (i.e., evolving with 
e -iwt^ m j xes creation a" and destruction operators. 

Using Eq. (21) and the normalization of Eq. (32) for 
the frequency eigenmodes V", it is easy to check that 
and a^} are creation and destruction operators, since 
their commutation relations, implied by the canonical 
commutators of the fields given in Eq. (7), are 

ala^] =i8{u-J)5w, (34) 

where now f3 can indifferently belong either to P or N. 

C. Matching conditions 

In this paper, we consider only configuration with a 
single analogue black hole horizons. For the sake of sim- 
plicity we model this system with two homogeneous half- 
line regions (representing, respectively, the interior and 
the exterior of the analogue black hole) connected by a 
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discontinuity at x — 0. On the two sides of the disconti- 
nuity, the oscillator strength takes different constant 
values 



A = /3i,LO(-x)+Pi,R0(x)- 



(35) 



We also impose that the inertia (/J,-f2f physically 
corresponding to the masses of the oscillator fields [see 
Eq. (1)], is the same in the two regions, i.e., f2i must 
varies accordingly to 



(36) 



In this optical system, the step-like profile has been 
proved to provide reliable result when the parameters 
do not vary very much between the two homogeneous 
regions [28] . This is indeed the case in those optical sys- 
tem, where the relative difference in the refractive index 
in the two regions is of the order of O.f % [9] (in the fol- 
lowing analysis, for illustrative purposes, we shall mainly 
use a much larger, but still quite small relativa difference 
of 10%). 

In this geometry, a frequency eigenmode V£ [see 
Eq. (20)] can be written as 

V° - E L " Kl K-x) + Kr ( 37 ) 

a ol 

where L" and are constant, and V" L and V" R are 
frequency-momentum eigenmodes, as in Eq. (22), that is 
they are solutions of the field equation Eq. (17) in the 
homogeneous left (x < 0) and right (x > 0) regions, 
respectively. 

The relations between L"'s and R"'s, are determined 
by solving the field equation in a neighborhood of x = 0. 
Writing V" of Eq. (37) in a more compact form 



V£ = V L 6(-x) + V R 0(x), 



(38) 



its first and second spatial derivatives are 

C = VI 9(-x) + V' R 9(x) + (Vr - V L )S(x) (39) 



C" = V'i 0(-x) + V' R ' 6(x) + 2(V R - V[)6(x) 
+ (V R -V L )5'{x). 



(40) 



We now put the above expressions in Eqs. (11), (12), (13), 
and (11), and group all the terms with 9(x), 9(—x), S(x), 
and S'(x). The terms with 9{x) and 0{— x) trivially give 
two sets of equations, separately satisfied by Vl and Vr in 
the homogeneous left and right regions. These equations 
are solved as in Sec. II B. For instance, Eq. (11) gives 



A R = 4ttc 2 U a ,r, A L = 4ttc 2 U a r. 



(41) 



The connection formulas between the two asymptotic 
regions are given instead by S(x) and S'(x) terms. Equa- 
tion (11) docs not contain any of them. Equation (12) 
has a term with S(x) 



v (Pi, R - P itL ) S(x) = 0, 



(42) 



yielding 

v[Pi,R(0)-P iiL (0)} = 0, (43) 

which has two solutions 

v = 0, or P ilR (0)=P ili (0). (44) 

Thus, when v ^ 0, Pi's must be continuous at x = 0. 

Equation (13) contains a second derivative of A, gen- 
erating both a term with S(x) and one with S'(x). The 
associated equations are 



A R (0) = A L (Q), A' r (0) = A' l (0), 



(45) 



so that both A and A' are continuous in x = 0, that is, the 
magnetic field B is continuous. Consequently, because of 
Eq. (41), also 11^ and 11^ are continuous in x = 0. Since 
the electric field is proportional to the time derivative of 
A, this implies that both E and its first spatial derivative 
are continuous in x = 0. 

Finally, Eq. (14) contains the first derivative of ITp (v 
is constant), generating a 5(x) term that yields 



v = o, or n Pi , L (o) = n Fi)fl (o), 



(46) 



and, consequently, ITp^s are continuous at x — when 
v ^ 0. 

Summarizing, we proved that 

• A, Ha, and their first and second derivatives are 
continuous in x = 0, 

• if v 5^ 0, Pi and ITp i are continuous. 

To completely fix the relations between modes in the left 
and right regions, one must determine the matching con- 
dition of P[ and H'p. , and also those of of Pi and 11^ if 
v = 0. 

a. u^O. The 9{— x) and 9{x) terms of Eq. (12) are 



"Pit = + P '- 



7 



' 'n Pi r-lA L 



7 Z V c 



(47) 



Since A, Ha, and P = iuP are continuous, and Pi^lf are 
equal on both sides of the discontinuity, Eq. (47) implies 
that also P- is continuous at x = 0. 

Similartly, the matching condition for Hp is found 
from the 6(— x) and 9(x) terms of Eq. (14): 



vW, 



Pi-L 



-Pi. 



fc,L (n 

fc,R (n Pi ,R - vW P . tR j = -p itR . 

By continuity of Pj 

A (n^-un'p.) = -p 



(48) 



(49) 



() 



is continuous at x — and, accordingly, lip. is discontin- 
uous. 

The continuity of P and P' implies that the derivative 
of P with respect to the laboratory time T 

drPi = ~f(Pi - vP!) (50) 

is continuous. Analogously, Eq. (49) implies that 

PidrlLpi (51) 

is continuous at x = 0. The physical meaning of this 
condition is evident when the second equation of (6) is 
rewritten as 



n 



7 



7 
c 



Differentiating with respect to T: 



d%P + ±d T A, 



Afi 



7. 

c 



(52) 



(53) 



where we used the constancy of ftfi 2 . Multiplying both 
sides of this equation by ft, using Eq. (49) in the form 



A 9 T n P = - 7 p it 



(54) 



and noting that P = —drA/c is the electric field, one 
obtains 

d\P = -f2?Pi + &fi 2 P. (55) 

This equation describes the dynamics of a harmonic oscil- 
lator of frequency fii and mass (/Sifi 2 )" 1 , which is forced 
by P. When the perturbation, caused by the laser pulse, 
passes at some X = Xq, at time T = Tq = Xq/v, the 
mass (Afi 2 )" 1 remains unchanged, but the oscillator fre- 
quency fli changes from fij iP to fi^z,- As a consequence 
of Eq. (55), P and <9 P P are unchanged as the perturba- 
tion arrives at X = Xq, but <9 T P has a finete jump, i.e., 
it is discontinuous at X — Xq and T = Tq. 

b. v = 0. In this case, Eqs. (12) and (14) simplify 

to 



(56) 



n 



p, 



ft 



Differentiating the first equation with respect to time, 
using the second equation and E = —A/c: 



Pi = -nf p + ftft? e, 



(57) 



which coincides with Eq. (55) because T = t when v = 0. 
Since PiVlfE is continuous at x = 0, 



P + nf p 



(58) 



is continuous and, for an cigenmodc of frequency uj 



oi r p, 



n? n) P 



(59) 



at x = 0. Thus, Pi is not continuous at x = 0, because 
f2i is discontinuous. 

Similarly, by derivating the second equation of (56) 
with respect to time, one obtains 



ri Pi = -n 2 t ( n Pi - 



which implies that 



is continuous, that is 



n 



Pi 



fi- 



ll 



Pi 



.L 



fi 2 



n 



Pi.R 



(60) 



(61) 



(62) 



at x = 0, and lip; is discontinuous. 

To conclude this section, it is worth checking how 
many constants are fixed by the matching conditions. In 
Eq. (37), there are 16 unknown constants L" and P", cor- 
responding, respectively, to the 8 solutions of the disper- 
sion relation (27) in each of the two homogeneous asymp- 
totic regions. Moreover, for every frequency cj, there are 
in general 8 global solutions V". Thus, those constants 
should be constrained by 8 equations, but we found 16 
matching conditions for the continuity of A, Pi, Ha, II Pi , 
and of their first derivatives. The system is apparently 
overdetermined. However, it is possible to show that 8 
constraints are redundant. In fact, the 8 continuity rela- 
tions of the conjugate momenta are directly implied by 
the 8 relations of the fields [see, for instance, Eq. (41)]. 

At the end of the day, for each frequency ui, there are 
16 parameters and 8 independent constraints, leaving 8 
free parameters, associated with the 8 globally defined 
solutions of the field equation (17). 



III. MODE ANALYSIS 

In this section we will apply the quantum field for- 
malism introduced in the previous section to build the 
global eigenmodes of the system by imposing the suitable 
matching conditions at the transition interface. Identi- 
fication of the different in-going and out-going channels 
and then calculation of the S matrix describing the scat- 
tering of light at the interface will be the crucial ingre- 
dients to calculate the intensity and the spectrum of the 
quantum vacuum emission in the next section. 



A. Asymptotic modes 

A stationary system (in the reference frame comov- 
ing with the laser pulse) made of two asymptotic ho- 
mogeneous regions, connected by a transition region was 
first investigated under a purely kinematic perspective in 
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Ref. [20]. In particular, it was shown that it is possible 
to tune the velocity of propagation of the pulse in such a 
way that the transition between the two asymptotic flat 
regions looks like an analogue black horizon. 

In that analysis, only the optical branch (correspond- 
ing to optical frequencies) of the dispersion relation was 
considered: within this approximation, the dispersion re- 
lation admits 4 solutions of k sharing the same value 
of the comoving frequencies lu. Before studying the full 
problem, taking into account all the branches of the dis- 
persion relation, it is worth summarizing the most rel- 
evant aspect of that analysis. In fact, the use of the 
Sellmeier dispersion relation introduces additional com- 
plications with respect to simpler subluminal dispersion 
relations usually considered in the literature [4]. In par- 
ticular, the absence of a low frequency regime forbid an 
immediate individuation of an analogue geometry. 

Assuming that the perturbation is moving in the pos- 
itive x direction (y > 0), a frequency dependent hori- 
zon is present for a giving comoving frequency lu when, 
on the left of the perturbation, there are only negative 
group velocity (measured in the comoving frame) modes, 
while on the right there are both negative and positive 
group velocity modes. In this case, light can propagate 
only leftward in the left region, both left- and right-ward 
in the right region. In analogy with black hole physics, 
the left and right regions correspond, respectively, to the 
interior and the exterior of a black hole, and, in the tran- 
sition region, there is one point corresponding to a black 
hole horizon. 

This physics is illustrated in Fig. 1, where the opti- 
cal branch of the dispersion relation (27) is plotted in 
the comoving frame, both in the left (left panel) and in 
the right (right panel) regions, for a step-like pulse (top 
panel) moving rightward at v — 0.66c. In the right region 
the dispersion relation is given directly by Eq. (27), while 
in the left region (representing the interior of a prop- 
agation pulse), the refractive index has been increased 
by Sn = 0.1. Solid (dashed) curves denote branches 
with positive (negative) laboratory frequency. As demon- 
strated in Sec. II B, they correspond to positive (nega- 
tive) norm modes. 

The dispersion relation is solved for a given comoving 
frequency lu. The arrows indicate the direction of prop- 
agation of the corresponding modes. Modes are named 
using the notation introduced in Sec. II. The superscript 
o stand for positive norm optical branch, o for negative 
norm optical branch. The subscript L and R denote, re- 
spectively, modes defined in the left and right region. For 
outer points (right panel, Sn = 0), the dispersion rela- 
tion has 4 solutions. Three of them correspond to modes 
(Kj R' K^.fl) propagating leftward from x = +oo 

toward the horizon. The fourth one, denoted by an empty 
dot, corresponds instead to an outgoing mode (Kj,_r), 
with positive group velocity. For inner points (left panel, 
Sn > 0.1) the dispersion relation has only 2 real solu- 
tions, both corresponding to leftgoing modes (V° L , V° L ) 
which propagate from the horizon toward x — — oo. Note 



that the two solutions present only in the right region for 
w m in < lu < Sai correspond, respectively, to a left- 
( v w%) and a right-going (V U<R ) mode. 

As a result, according to the above given definition, 
within the frequency range 0J m i n < lu < w max (where 
tkViaia and Sw max are represented in the right panel with 
dotted lines), the discontinuity in the pulse profile rep- 
resent an analogue black hole horizon. Nota that there 
is one negative norm mode V° R , leftward propagating 
from +oo to the horizon, and one solution (right panel, 
empty dot), corresponding to a mode V u ,r propagating 
from the horizon to infinity in the exterior part of the 
analogue black hole. This mode structure might origi- 
nate some phenomenon similar to Hawking radiation. 

In the present paper, this analysis is extended to the 
full Sellmeier dispersion relation, which in general admits 
8 solutions, corresponding to 8 (propagating if k is real) 
asymptotic modes (AM) in each region (x < or x > 0) . 
Proceeding in analogy with Ref. [20] , one must first iden- 
tify a range of frequencies Lu m - ln < lu < w max , for which 
the optical branch possesses 4 real solutions in the region 
and 2 in the left one. Second, on each side one must lo- 
calize 4 additional solutions on the other branches of the 
dispersion relation and describe the associated modes. 

In Fig. 2 the full Sellmeier dispersion relation is plotted 
for x < (upper left panel) and x > (upper right panel) 
in the laboratory reference frame (fl,K) [see Eq. (29)]. 
A boost is then performed on the axes and the new axes 
(lu, k) in the reference frame comoving with the pulse are 
drawn [see Eq. (27)]. The central region of those plots 
(gray dot-bordered square) is zoomed in the bottom pan- 
els. The dispersion relation is graphically solved for a 
fixed value of the comoving frequency lu (dashed line), 
with w m i n < lu < Lu max . There are 8 branches: 4 with 
positive (solid curves) and 4 with negative fi (dashed 
curves) , symmetrically placed in the lower plane. In this 
configuration lu is small enough that no solution belongs 
to the highest (positive or negative) energy branches. 
We therefore name only the 6 branches of the disper- 
sion relation with low energy |f2|. In the upper plane 
(f2 > 0), starting from the lowest energy branch we call 
them lower (1), optical (o) and upper (u). Symmetri- 
cally, the 3 branches with negative laboratory frequency 
n and negative norm (as demonstrated in Sec. II C) are 
labeled by 1, 6, and u. Accordingly, the solutions of the 
dispersion relation arc labeled by a superscript 1, o, u, 
1, 6, and u. To ease the reader the comparison with 
Fig. 1, the solutions on the positive and negative fre- 
quency branches o and 6 are denoted by an empty dot. 
The arrow above each solution indicates the direction of 
propagation (group velocity in the comoving frame) of 

the associated mode V"^ R . 

Note that, for x < (left panels), there are 6 real- 
k solutions, all corresponding to leftgoing modes. The 
remaining two solutions of Eq. (27) have complex conju- 
gate k. They are associated with exponentially growing 
(V^f£ w ) and decaying (V^ c £) modes for x — > — oo. For 
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FIG. 1. Simplified analysis of the Sellmeier dispersion relation (27) in fused silica, as seen from the comoving frame. Only 
the optical branch is shown. Positive (negative) laboratory frequency branches are represented by solid (dashed) curves. The 
dispersion is plotted on the left (Sn = 0.001, left panel) and on the right of a perturbation (Sn = 0, right panel) moving 
rightward with V = 0.66c in the laboratory frame (see top panel). The black horizontal line represents a generic frequency for 
which there are 4 real solutions in the right region (right panel), 2 real solutions in the left one (left panel), and the system 
shows an analogue black hole horizon. In the right panel, the dashed horizontal lines indicate the maximum and the minimum 
values of the frequency for which this behavior occurs. Modes are labeled with the notation introduced in Sec. II and their 
propagation direction is indicated by arrows. 



x > (right bottom panel), instead, the 8 solutions are 
all real. The two extra real solutions, that do not have a 
corresponding solution on the left side, belong to the op- 
tical branch and are associated, respectively, with a left- 
going mode, named V° 2 R , and with the unique rightgoing 
mode, simply named V u< r, without any superscript. 



B. Globally defined modes 

In the previous section, plane waves modes propagat- 
ing in the asymptotically flat left and right regions have 
been identified. Combining those asymptotic waves, two 
relevant bases of globally defined asymptotically bounded 
modes (GDMs) (not diverging at infinity) can be con- 
structed. 

We define the in basis as the set of in modes, whose 
asymptotic decomposition (37) has only one AM with 
group velocity v g directed toward x = 0. We say that the 
group velocity of an AM is directed toward the horizon 
if v g > (v g < 0) for modes which are solution of the 
mode equation in the left (right) region. 

Analogously, we define the out basis as the set of out 
modes, whose asymptotic decomposition have only one 
AM with group velocity directed toward x — — oo (x = 
+00), if the AM is a solution of the field equation in the 
left (right) region. 

To clarify the meaning of these definitions with an ex- 
ample, the asymptotic decomposition of one of the three 
negative norm in GDMs is schematically represented in 
Fig. 3. We name it V^ n, °, since its unique incoming AM is 
V° R , that is i?° is the only non- vanishing coefficient as- 



sociated with a mode with group velocity directed toward 
the horizon, in the asymptotic decomposition of Eq. (37) 
(see also Fig. 2). 
In Table I, all the coefficients L a 's and Rf's are re- 

' UJ UJ 

ported for the asymptotic decompositions of the three 
negative norm m GDMs V^ 5 , V^. As noticed 

in the previous section, both the AMs V°\ and V w r, 
which are present only in the right region, lie on the op- 
tical branch of the dispersion relation. As mention at 
the beginning of Sec. Ill A, in the simplified analysis re- 
stricted to the optical branch, analogue Hawking radia- 
tion is expected on the AM Vu,r, due to the scattering at 
the analogue horizon of the optical negative-norm incom- 
ing GDM V^ n, °. Thus, aiming to extend the definition 
of Hawking radiation to the present situation, adopting 
standard notation, we name /3 U the coefficient of V u> r in 
the expansion of V^ n, °. Similarly, we name a u the co- 
efficient of V° L . Furthermore, adopting the notation of 
Ref. [29], the coefficients of modes with negative (posi- 
tive) norm, in the expansion of modes with negative norm 
are named A (B). As an example, in the expansion of the 
negative norm mode V™'°, the coefficients of the nega- 
tive norm modes Vr, T and V" T are named A ; 1 and ^4°, ,u . 

UJ ,Li UJ ,/j UJ UJ 

The coefficients of the positive norm modes V° L , L , 
and V£ L are named B u , B°> 1 , and B°' u , respectively. 

Finally, the 8 unknown coefficients in the two left and 
right asymptotic expansions for each GDM are deter- 
mined by imposing the matching conditions (8 indepen- 
dent equations), derived in Sec. II C. 

By proceeding along this line, it is possible to con- 
struct 7 in GDMs: V£ n '\ V™> 5 , V™'", V£ n >°, V^ n < u , 
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FIG. 2. Graphical representation of the Sellmeier dispersion relation, as seen from the laboratory reference frame (il,K), for 
x < (left panels) and x > (right panels). The (oj, k) axes of the comoving reference frame are obtained through a boost of 
velocity v. As sketched in the top panel, the refractive index in the left region is larger than in the right region. This difference 
in the refractive index is obtained by properly changing the parameters Pi and Qi appearing in the Lagrangian (1). In this plot 
the values of the velocity v and of the refractive index change 5n have been arbitrarily chosen for illustrative purposes. The 
bottom panels are zoom of the gray dot-bordered squared of the respective upper panels. The dispersion relation is graphically 
solved for a fixed comoving frequency U), chosen in the frequency window in which the black hole horizon is present. Solutions 
appear both on the positive-norm positive-fi branches (solid curves) and on the negative-norm negative-Si branches (dashed 
curves). The empty dots denote solutions on the optical branches with positive (o) and negative (5) frequency £1. The arrows 
indicate the direction of propagation (group velocity in the comoving frame) of the associated modes V"^ R . In the left region 
(left panels), the dispersion relation has only 6 real-fc solutions. In the right region (right panels), the real-fc solutions are 8, 
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TABLE I. Coefficients of the asymptotic decomposition of the globally defined in modes Vj 1 ' , VZ n '°, and VZ n ' u on the bases of 
left and right AMs. 




of GDMs. 



V 



dw e 



J2 Vj?' a a%> a + ]T Vj?> & a^ 



\a£P 



a£N 



h.C. 



(63) 




^ o u t . a u t , a I 



h.C. 



(64) 



The transformation between the two bases follows 
straightforwardly from the construction of the previous 
section. In fact, since the two bases of in and out GDMs, 
the basis of left AMs and the basis of right AMs are all 
orthogonal, the coefficients of the matrix connecting the 
in and out bases can be directly read from the expansion 
of in and out modes on AMs. For instance, from the 
second line of Table I: 



V° - 



FIG. 3. Asymptotic decomposition of the globally defined in 
mode VZ n '°. The only asymptotic branch with group velocity 
directed toward the horizon is V£, R . Positive (negative) norm 
modes are represented by solid (dashed) lines. 



V^ n '° 2 , whose unique branches with group velocity di- 
rected toward the horizon are, respectively, V* R , V° R , 
V» R , V^ R1 VZ tR ,V» R , V°%. Similarly, one constructs 7 

out GDMs: V° ut ~\ V° ut > 5 , V^ out ' a , V° ut '\ V^ out '°, V° ut ' u , 
V° ut , whose unique branches with group velocity directed 
to infinity are, respectively, V^ L , V° L , V* L , V^ L , V° L , 

v U) r. 



C. The scattering matrix 

As in Eq. (33), the field operator V is expanded indif- 
ferently with respect either to the in or to the out basis 



+ B°'V 



It rout, 1 

UJ ^ UJ 



yyin.o . /io,lT^out,l _j_ Trout, 6 _j_ yio.UT7-out,u 

V uj -^uj uj ' ^uj " UJ T" ^ w *oj 

out.o i DO,UT/-out,u I a t rout 



+ B o,u y out,u + ^ F out_ (65) 



Repeating this procedure for each in mode, the scattering 
matrix S is fully determined 



(66) 



where /3 and fj' run over all positive and negative norm 
modes. The relation between in and out destruction and 
creation operators is easily derived from the S matrix: 



A out = S T A in 

where A ln 

i in = (^ lJt fifrM a£' St a; 



(67) 



in.l Ain,o Ain,u ~ in,o2 

Uj,. Uj. . Uj. . 



(68) 



£out — / "out, If ^out.of "out.uf "OUt.l Aout.o ~out,u "out 

A — \ a u ' a u ' a U) ' a u a u a uj a LO 

(69) 

are 7-dimcnsional vector formed, respectively, by the in 
and out creation operators a™' 01 ' and a™' 01 ' of the neg- 
ative norm modes, and by the in and out destruction 
operators a}£ ,a and a° ut: " of the positive norm modes. 
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IV. SPONTANEOUS EMISSION 

In this section we shall make use of the S matrix cal- 
culated in the previous section to provide quantitative 
predictions for observable quantities such as the inten- 
sity and the spectrum of the quantum vacuum emission. 
A straightforward extension of the formalism following 
the lines of [30, 31] can be used to obtain information on 
the correlation properties of the emitted quantum vac- 
uum emission into the different out-going modes: given 
the importance of correlation experiments to assess the 
quantum vacuum nature of the emission, this problem 
will be the subject of future work. 



A. Comoving frame 

For the sake of simplificy, we make the reasonable as- 
sumption that there are no ingoing particles, i.e., the 
system is in the Fock vacuum state defined by the de- 
struction operators associated with the in modes: 



lOi,, 



0, 



10, 



0. 



(70) 



The occupation numbers of out modes on this state are 
easily computed using Eq. (67). They in general do not 
vanish, since Eq. (67) mixes creation and destruction op- 
erators. For instance, the expected occupation number 
of the unique rightgoing mode V° ut on the state |0j n ) is 



(0i 



|27ra° utt a° ut 



|0 in )=2^(0)(|^| 2 + |^| 2 + |^| 2 ), 



( 71 ) 

where the factor 2-7T has been inserted coherently with 
the normalization of the Fock operators of Eq. (34) . As 
usual, there is an infrared divergence associated with the 
quantization of a field theory in an infinite spacetime 
volume. Going to a finite size time-box At, one must 
replace 



2tt(5(w - J) 



At 5^ 



(72) 



This implies that the number of particles An created in 
a time At at a frequency ui is 



An" 



(73) 



Analogously, there are several other channels in which 
particles are created, related to the mixing of the other 
positive and negative norm modes. For the positive fre- 
quency modes one obtains 



ArC U = (|£>| 2 
An° ut '° = (|£U 2 ■ 
ArC" u = {\Bl\ 2 +\B 



\B°J\ 2 
\B l f\ 2 

1,U|2 



\B n J\ 2 )At, 
\B n :°\ 2 )At, 



\B 



o,ui2 



)At. 



(74) 
(75) 
(76) 



Since in the comoving reference frame, the source of 
photons (the pulse) is at rest, the number of created par- 
ticles r w per unit time and unit bandwidth coincides with 



the flux of particles I u crossing a certain surface at con- 
stant x. Thus, the flux of particles per unit time and unit 
bandwidth in the comoving reference frame is 



r 



dn° 



A? 



dt duj 



At 



(77) 



The flux of energy (which coincides with the energy pro- 
duction rate in the reference frame where the source is 
at rest) associated with the mode a is 



J pout, a 

Ju dtduj 



(78) 



In Fig. 4 we plot the fluxes of particles I° ut (solid line), 
I™*' 1 (dashed line), 7° ut >° (dotted line), 7° ut < u (dot- 
dashed line) and the respective energy fluxes (right panel) 
fZ n \ fZ nt, \ fZ ut '°, and /° ut ' u , for a pulse moving with ve- 
locity v — 0.66c. The values of and f2j R in the right 
region have been chosen accordingly to the dispersion re- 
lation in fused silica, the material used in the experiment 
of Ref. [9]: 



P hR = 0.07141914, 
P 2 ,r = 0.03246304, 
fo R = 0.05539915, 



m ltR = 0.125285 eV, 
m 2 , R = 10.6661 eV, (79) 
m 3R = 18.1252 eV. 



In the left region, for illustrative purposes we used 

Pi, L - (i + e)ft >fl , n hL = (i + e)-v 2 n iiR , (80) 

with a very large value of e = 0.3. When f2 is in the opti- 
cal range and far enough from the poles of the dispersion 
relation 



Stir 



1 



n R 



2n R 



(81) 



that yields a quite large value of Sn « 0.12, which re- 
quires a very strong laser intensity I w 3 x 10 14 W/cm 2 [9, 
20]. 

From Fig. 4, the dominant contribution to the energy 
flux comes from the rightgoing modes, corresponding to 
a Hawking-like process, in the sense discussed in the pre- 
vious section and in Ref. [20]. What can we say about 
its thermal properties? It is well known that, the ther- 
mal energy flux in 1+1 dimension is exactly constant in 
the limit U) — > 0. In the present situation, however, this 
channel is open only in a narrow range of frequencies 
w m i„ < uj < uj max (see Fig. 1), and the low frequency 
limit cannot be taken. Therefore, asking whether the 
spectrum is thermal or not is not really a well-posed ques- 
tion. Nevertheless, one may observe that the energy flux 
is almost constant, far from the boundary values w m in 
and w max . This is all what is possible to say about the 
thermal properties of this spectrum. 

To investigate more deeply the features of the spectrum 
of V^ ut (the global modes associated with the asymp- 
totic plane wave V UtR , see Fig. 2), in Fig. 5 the flux of 
particles is separated into their components due to |/3o>| 2 
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FIG. 4. Particle fluxes (left panel) 7° ut (solid line), 7° ut '' (dashed line), J° ut '° (dotted line), J° ut ' u (dot-dashed line) and the 
corresponding energy fluxes (right panel) /° u \ /° ut '\ /° ut '°, and /° ut ' u as seen from the comoving reference frame, for a horizon 
configuration, as sketched in Fig. 2. Parameters v = 0.66c and e = 0.3, yielding via Eq. (80) a refractive index jump of Sn = 0.12 
for optical frequencies. 



(solid line), \/3l\ 2 (dashed line) and \/3^\ 2 (dotted line). 
Note that the contributions to spontaneous emission on 
the mode V° ut by the lower and the upper branch are 
negligible with respect to the term coming from the op- 
tical negative frequency mode V^ n, °, associated with the 
asymptotic mode V° R of Fig. 2. This behavior was pre- 
dicted and justified in Ref. [20] and it has been now con- 
firmed by the present calculation. 



B. Laboratory frame 

One may wonder what are the spectral features as mea- 
sured in the laboratory reference frame. Note that, the 
relevant experimental quantity is not the flux of particles 
but the total number of particles created by the moving 
pulse. In fact, the pulse stably propagates only for a very 
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FIG. 5. Separate contribution of |,(3^| 2 (solid line), l/^j 2 
(dashed line) and |/3^| 2 (dotted line) to the flux of particles 
I^ ut , for the horizon configuration considered in Fig. 4. 



short distance AX S (about 1 mm [10]), corresponding to 
a duration of AT S — AX s /V g . In an experiment one 
observes the pulse in this small region and measures the 
photons produced therein. Naming Rq the production 
rate of particles at frequency fi, as measured in the labo- 
ratory, the number of particles produced by the pulse 
in the laboratory time dT, in the range of frequency 
(Q,Q + dfi) is 



AN = R n dTdn. 



(82) 



Since the number of created particles is invariant under 
Lorentz transformation 



AiV = r n dtduj. 



(83) 



that is AiV can also be computed as the number of par- 
ticle created in the comoving time dt by the pulse, in the 
frequency range (uj,u> + duj). To compute the relation 
between Rq and r w , it is enough to compute the trans- 
formation rules of the frequency range and of the time 
interval. Using the inverse of the Lorentz transforma- 
tion (28) 



lu = 7(f2 — vK) 
du) can be expressed in term of d£l as 



dK 

du> = 7 ( dfl — v——-dfl 
dil 



7 1 



V g {Q) 



(84) 



dil, (85) 



where V g is the particle group velocity measured in the 
laboratory reference frame. The tranformation of the 
time interval is found by noting that the source (the 
pulse) is at rest in the comoving frame, so dt is its proper 
time interval. Consequently the corresponding labora- 
tory time interval dT is given by the usual Lorentz dila- 
tion of time 



dT = jdt. 



(86) 
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FIG. 6. Energy production rate £Q Ut on the rightgoing mode 
V°(~Qj as seen by an observer in the laboratory, for the horizon 
configuration considered in Fig. 4. 



Putting everything together, the rate of particle produc- 
tion, as seen in the laboratory frame is 

Rfi =(l-~)r u . (87) 

The corresponding energy production rate en in the 
laboratory frame is 

e n = nrtR Q = mi a (i - ~) ■ ( 88 ) 

This quantity is plotted in Fig. 6 for the rightgoing 
mode V° ut . The effect of the boost is twofold. First, 
it shifts the emission frequency in the laboratory: £1 is 
much larger than the comoving frequency oj, reaching op- 
tical frequencies. Unfortunately, as observed in [20], the 
boost effect is lost when one looks perpendicularly to the 
pulse propagation in a 3-dimensional system. Radiation 
at such high frequency could be observed perpendicularly 
only if some scattering process changed the direction of 
photons. 

To conclude this section, let us roughly estimates the 
number of photons emitted on the rightgoing mode 
by this process. Integrating in frequency the production 
rate in the laboratory frame and multiplying it by the 
time AT S = AX S /V over which the pulse stably propa- 
gates, we obtain (for AX S ps 1 mm) 

/AY 
RncM ps 240. (89) 

Of course, the same result can be obtained by directly 
integrating the particle production rate r u , measured in 
the comoving frame (coinciding with the flux I u ), over 
the comoving frequency oj and multiplying by the proper 
duration of the pulse At s = AT s /j(V). 



C. Outside the analog horizon frequency window 

In this section, the emission properties are investi- 
gated for comoving frequencies below w m i n and above 
w max , where the main particle production channel dis- 
appears. Particle creation is indeed possible also outside 
the Wmin -7- oj max range, through the other channels and 
the dashed, dotted and dot-dashed lines of Fig. 4 can 
then be extended for uj < oj m i n and oj > oj max . Since, for 
those values of the comoving frequency, there are only 
6 real-fc solutions also for x > 0, a new mode analysis 
must be performed following Sec. Ill to compute a new 
6x6 scattering matrix. In Fig. 7 the numerical results of 
this computation are reported, together with the results 
for w m i n < oj < oj max , obtained in the previous section. 
The flux of particles in the comoving frame /£ ut '' (dashed 
line), 7° ut '° (dotted line), and I° ut < u (dot-dashed line) are 
plotted in the left panel, and the respective energy fluxes 
fZ n \ fZ nt '\ fZ ut '°, in the right one. 

First, note that the emission is peaked in the frequency 
range cj m j n < oj < w max , thanks to the presence of the 
Hawking channel (solid line). This result confirms the 
naive expectation that the presence of horizon should 
enhance the production of particles. 

Second, both the occupation numbers and the energy 
fluxes are continuous at uj = u) m in and ui = oj maK , for 
modes that do not feel the presence of any horizon, i.e., 
when no turning point is present. For these modes the 
transition between the two regimes is continuous. 

Third, it is worth remarking that the high frequency re- 
gion of this plot must be read cum grano salts. Indeed, we 
used step functions to describe the spatial behavior of 
and Oj, such that arbitrary large frequency/momentum 
modes are excited. However, in real physical situations, 
the transition between the two regions x < and x > 
takes place on a finite length and only modes up to a cer- 
tain frequency w cu toff are excited. Thus, Fig. 7 provides 
reliable results only up to w C utoff- Beyond this frequency, 
the off-diagonal coefficients of the scattering matrix 5 go 
to zero exponentially with ui and S reduces to the iden- 
tity. Indeed, for uj > a; cuto ff, modes do not mix because 
they are well approximated by their WKB expansion. In 
a realistic situation, w C utoff is not larger than the fre- 
quency associated with the width of the steepening of 
the pulse [10, 16] (approximately 1 (an), corresponding 
to ^cutoff ps 2 x 10 15 s {hu cutoS ps 1.3 eV). 



V. PARTICLE PRODUCTION IN OTHER 
CONFIGURATIONS 

In the previous section, we derived the flux of spon- 
taneously created particles in a configuration closely re- 
sembling a black hole geometry. However, radiation from 
vacuum fluctuation is expected even in optical system 
without horizons, provided that negative norm modes 
are present [20]. We now apply the techniques intro- 
duced in this paper to the horizonless configuration ex- 
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FIG. 7. Fluxes of particles 7° ut (solid line), /° ut '' (dashed line), 7° ut '° (dotted line), J° ut ' u (dot-dashed line) and the respective 
energy fluxes (right panel) /^ ut , /<£ , /^ ut '°, and fZ ut ' u , for the horizon configuration considered in Fig. 4. The plot has been 
cutoff at cj > Wcutoff- IZ ut an d f^ ut (solid lines) are defined only for w m in < cj < w m ax. 



perimentally realized in Ref. [9] for a weaker perturbation 
(Sn ps 0.001). Indeed, outside that frequency range, the 
two optical modes V^ ut and V^ ut '° 2 (see Fig. 2) disap- 
pear. Thus, this situation is completely analogous to a 
system without horizons. As a preliminary step, however, 
we compute the emitted flux for a horizon configuration 
with a value of Sn ~ 0.001, which coincide with the ex- 
perimental value [9], so that it is possible to compare 
the flux produced by a weak perturbation, both in the 
presence and in the absence of horizons. 



A. Small refractive index jump 

In a realistic situation Sn is generally a couple of or- 
ders of magnitude smaller than the value above used. 
Unfortunately, if Sn is small, the pulse velocity v must 
be extremely fine tuned in order to obtain a configura- 
tion with a horizon, as in Fig. 2. The frequency win- 
dow (u min ,w max ) becomes very narrow and the emission 
on the rightgoing mode V^ ut is strongly suppressed. In 
Fig. 8, the energy fluxes associated with the spontaneous 
particle production are represented for Sn = 0.001. This 
value has been obtained by modifing the values of (3i and 
Qi in the left region, as in Eq. (80), with e w 0.0026. The 
pulse velocity is v — 0.6838c, hoj m - ln = 0.01379 eV, and 
hw max = 0.01396 eV. 

The Hawking channel is not the dominant one, but its 
flux is now comparable with the flux of leftgoing particles 
on the optical-branch mode V^ ut ' . In Fig. 9, the fluxes 
of particles and the respective energy fluxes are plotted 
for any frequency smaller than the cutoff w cut off, arbi- 
trarily chosen at huj cnto s ~ 1.3 eV. Note that the range 
(w m i n , w max ) is so narrow to be almost invisible. As a 
consequence its contribution to the total flux of leftgoing 
particles is quite small, although the maximum intensity 
is reached in this range. 

Furthermore, the number of particles produced on the 
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FIG. 8. Energy fluxes /° ut (solid line), /° uU (dashed line), 
/° ut '° (dotted line), and /° ut ' u (dot-dashed line), for a con- 
figuration with a small refractive index jump (Sn w 0.001), 
obtained with e = 0.0026 [see Eq. (80)]. The pulse velocity 
has been fine tuned to v = 0.6838c, so to have an analogue 
horizon. 



lower and the upper branches are about 3 order of mag- 
nitude smaller the the production on the optical branch 
(both right- and left-going). We have also checked that 
the main contribution to particle production comes from 
the incoming negative frequency mode on the optical 
branch. Both those results confirm that the lower and 
upper branch can be safely neglected as anticipated in 
Ref. [20]. Note that the latter result was shown in Fig. 5, 
for a different configuration with larger refractive index 
jump. 

Finally, in Fig. 10, the energy produced per unit band- 
width and unit time on the rightgoing branch V£j&\ is 
plotted in the laboratory reference frame. Quite surpris- 
ingly the frequency windows is much larger. In fact, even 
if w m in and w max almost coincide, k varies in a wider 



15 



Ha)- 



fia> c 



0.001 



10- 



10" 



10" 



10" 




0.001 



0.01 0.1 



10" 3 
_ 10" 6 
"a 10" 7 
10" 8 
10" 9 







- 
- 




' / 






. / 






- ^ - . - ' " / * -> 






■ / 











0.001 



0.01 0.1 
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FIG. 10. Energy production rate £n ut of the rightgoing mode 
V^m) as seen by an observer in the laboratory, for the horizon 
configuration considered in Fig. 8. 



range, because the group velocity v g , measured in the 
comoving frame, almost vanishes. Consequently, from 
Eq. (28), fl varies in a quite wide window. However, 
since the total number of created particles must be the 
same in both reference frames, the production rate per 
unit bandwidth is much smaller. Integrating the pro- 
duction rate over frequency and multiplying for AX S /V, 
with AX S « 1 mm, we obtain a tiny average number 
(N « 8 x 10 -4 ) of particles produced by each pulse. In 
this configuration, fewer particles are produced on the 
rightgoing mode V° ut , with respect to other modes. In 
fact, the total number of photons produced on positive 
norm modes, obtained by integrating all the curves of 
Fig. 9, left panel, is N tot « 7 x 1(T 2 . 



B. Horizonless configuration 

A system without any horizon is similar to the case 
investigated in Sec. IV C for uj < Lo m - m or uj > w max , with 
the exception that there is now no frequency window with 
8 real-fc solutions of the dispersion relation in the right 
region. This configuration can be obtained, as in Ref. [9], 
by making the pulse so fast that the two rightest solutions 
on the optical branch disappear (see Fig. 1, right panel, or 
Fig. 2, bottom right panel). The pulse velocity is tuned 
so that the straight dotted line of the bottom panel of 
Fig. 2 becomes steeper than the tangent to the optical 
branch at its inflection point. 

In analogy with Fig. 1, in Fig. 11 the dispersion rela- 
tion in the comoving frame is solved only for the optical 
branch, on both sides of the perturbation, for a pulse with 
Sn = 0.001 and moving at v w 0.69c, as in the experi- 
ment of Ref. [9] . Only two solutions (one with positive 
and one with negative norm) are present both in the left 
(V° tL , V° >L ) and in the right V% R ) regions for all 

values of the frequency w. In particular, the solutions 
corresponding to AMs V° 2 R and V^^r, the mode respon- 
sible for analogue Hawking radiation disappear. Even if 
particle production is still possible, due the the mixing of 
positive and negative norm frequency modes, the spec- 
trum will be different and the Hawking-like channel is 
absent. For the sake co completeness, in Fig. 12, the full 
Sellmeier dispersion relation is graphically solved for this 
configuration, for x < (left panels) and x > (right 
panels) . 

Finally, the fluxes of emitted photons (Fig. 13) are 
suppressed by the smallness of Sn, as in Sec. VA. The 
average number of created particles for a pulse which 
stably propagates for about 1 mm is indeed ps 0.06. 
Moreover, the only effect caused by the absence of the 
horizon is the disappearance of the outgoing rightgoing 
mode. The integrated (over frequency) number of cre- 
ated particles does not significantly differ with respect to 
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FIG. 11. Simplified analysis of the Sellmeier dispersion relation (27) in fused silica, as seen from the comoving frame. Only 
the optical branch is shown. Positive (negative) laboratory frequency branches are represented by solid (dashed) curves. The 
dispersion is plotted on the left (Sn = 0.001, left panel) and on the right of a perturbation (Sn = 0, right panel) moving 
rightward with V = 0.69c in the laboratory frame (see top panel). For all frequenies there are 2 real solutions both in the 
left (left panel) and in the right region (right panel). Modes are labeled with the notation introduced in Sec. II and their 
propagation direction is indicated by arrows. 



the horizon configuration with the same small vlues of 
Sn, because of the narrowness of the frequency window 

k^min ~ ^max- 

For the sake of simplicity, the perturbation has been 
described by a step-like profile. The propagation of light 
propagating in such optical medium has been modeled 
by the Lagrangian of the electromagnetic field coupled 
with polarization fields, in a one-dimensional system. 



VI. SUMMARY AND DISCUSSION 

In this paper we have developed a microscopic theory 
of the spontaneous quantum vacuum emission generated 
by a strong light pulse propagating in a Kerr nonlinear 
optical medium: the effect of the pulse is modeled as a 
moving refractive index perturbation following in a lo- 
cal and instantaneous way the pulse. Here we restrict 
our attention to the case of a single sharp interface sep- 
arating two homogeneous regions of spatially constant 
optical properties for which several authors have antici- 
pated the occurrence of the optical analogue of a black 
hole horizon for suitable values of the pulse speed and the 
amplitude of the refractive index jump. With respect to 
previous work, our theory takes into full account the un- 
avoidable frequency refractive index dispersion of the di- 
electric medium and identifies the different regimes with 
or without horizon that can be obtained depending on 
the pulse parameters. 

Moving to the reference frame comoving with the pulse 
where the optical properties of the system are time- 
independent, the classical eigenmodes of the coupled elec- 
tromagnetic and matter polarization fields at a given fre- 
quency are derived within a Lagrangian formalism and 



then quantized through canonical quantization. The 
mixing of positive and negative norm modes at the same 
frequency is responsible for the emission of quantum vac- 
uum radiation, whose rate of production has been com- 
puted. 

In a configuration with horizon, one can identify a 
quantum vacuum emission channel showing some re- 
markable similarities with Hawking radiation: this emis- 
sion channel completely disappears in the absence of hori- 
zon and quasi-thermality of its emission spectrum is sig- 
naled by an almost constant flux of energy emission when 
observed in the comoving frame (of course, this spectrum 
turns out to be distorted when observed in the labora- 
tory frame). On the other hand, because of the nontrivial 
shape of the Sellmeier dispersion relation, this emission 
channel is only available and active for frequencies in the 
restricted uj m - m < lo < cj max range where the analog black 
hole horizon exists. In addition to this Hawking-like ra- 
diation, we have found that quantum vacuum radiation 
is also emitted on several other channels and outside this 
frequency window, albeit with a much weaker intensity, 
a strongly non-thermal spectral distribution, and inde- 
pendently of the presence or not of the horizon. Future 
work will extend our theory to regimes showing analogue 
white hole horizons. 

As the experiment in [9] was carried out in a param- 
eter range not showing any horizon, it is interesting to 
conclude the paper by discussing whether or not it is 
legitimate to denote the quantum vacuum emission in 
this regime as "Hawking radiation" . This question is all 
more relevant given the on-going debate [14, 15] on the 
interpretation of the experimental results. Even if the 
results of our calculations are objective, answer to this 
question requires a preliminary agreement on the defini- 
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FIG. 12. Graphical representation of the Sellmeier dispersion relation, as seen from the laboratory reference frame (f2, K), for 
x < (left panels) and x > (right panels). The chosen plot range excludes the highest frequency branch. However, there are 
no solution on it, see discussion in Sec. Ill A. The (to, k) axes of the comoving reference frame are obtained through a boost of 
velocity v. As sketched in the top panel, the refractive index in the left region is larger than in the right region. This difference 
in the refractive index is obtained by properly changing the parameters Pi and Qi appearing in the Lagrangian (1). In this plot 
the values of the velocity v and of the refractive index change Sn have been arbitrarily chosen for illustrative purposes. The 
dispersion relation is graphically solved for a fixed comoving frequency cj. Solutions appear both on the positive-norm positive- 
Si branches (solid curves) and on the negative-norm negative- branches (dashed curves). The empty dots denote solutions on 
the optical branches with positive (o) and negative (5) frequency Q. The arrows indicate the direction of propagation (group 
velocity in the comoving frame) of the associated modes V^^ R . The dispersion relation has only 6 real-fc solutions in both 
regions. 



tion of Hawking radiation, which is somehow a matter of 
personal taste. 

If a broad definition is chosen, where the only require- 
ment is the steady production of particles out of the 
quantum vacuum, there are no difficulties in consider- 
ing the predicted emission as an example Hawking radi- 
ation. The answer is different if we also require the emis- 
sion to be thermal, which is definitely not the case of the 
quantum vacuum emission in horizonless configurations 
or outside the w m i n < ui < w max window of a horizon 
configurations. And even inside the uj m - ln < u < w max of 
a horizon configuration some care has to be paid not to 
overstretch the gravitational analogy: even if the emis- 
sion has an approximately thermal spectrum within this 
window, the fact that is does not extend down to low 
frequencies is related to the difficulties of a description 
of the light propagation in terms of a curved space-time 
metric as done in other analogue models. 
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Appendix A: Lagrangian density in the comoving 
frame 

In Sec. II A, we applied a Lorentz boost A to the sys- 
tem described by the Lagrangian of Eq. (1), to obtain 
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Eq. (4). In doing so, we transformed the space and time 
coordinates, but we treated both the fields A and P as 
scalars. 1 This is possible since we are dealing with a 1+1 
dimensional system. 

To properly proceed, one should first apply the Lorentz 
transformation A to the fields A and P, assuming that 
both A and P oscillate in the z direction. A is the z 
component of the electromagnetic potential 



A^ 



which transform as 



o 



W 



(Al) 



(A2) 



where the subscript p indicates that A£ is measured 
in the reference frame comoving with the pulse. Anal- 
ogously, P is one of the three components of the 
magnetization-polarization tensor 
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-cP x 

-CPy 

\-cP z -M, 



cP x 
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-M z 


M z 
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My 

-M x 
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which transform as 



p 



M pa . 



( cP\ 


\-cP / 
(A3) 

(A4) 



By doing so, the Lagrangian density can be rewritten in 
term of the new fields Ag and Mg v . 



1 To maintain the notation compact, the subscript i is omitted in 
this Appendix. 



Moreover, the Lagrangian density can always be ex- 
pressed by using new fields that are functions of A^ and 
M£ v , provided that the transformation to those new 
fields does not involve any time derivative of Ai^ and 
Mg v . Such a transformation generates a canonical trans- 
formation on the Hamiltonian variables, so that the com- 
mutation rules among the new fields and their respective 
momenta are still canonical. We chose to apply to A£ 
and M!p V the linear transformation A -1 , not involving 
any time derivative of the fields. The new fields A£ and 
M£" are 

A p = (A" 1 )^ = (A-y„A V = A p 
M^ = (A- l r p (A- 1 r a MP' (A5) 
= {K- l Y p {K- 1 Y a A p a K' 1 pM^ t3 = AP". 

That is, A% and A/#" coincide with the old fields A^ and 
AP" ', as measured in the laboratory frame. The only 
nonvanishing components of those fields are A and P, as 
in Eqs. (Al) and (A3). 



Appendix B: Mode normalization 

The scalar product (18) on two eigenmodes with real 
frequencies uj\ and W2, and real momenta k ai and k a2 is 



(V ai V a2 \ = - I rlT p -i^2-^l)t+i{k a2 -k ai )x T/ait „ T/<*2 
\ v lu 1 i v cj 2 I hi 1 2 

= 2n6(k a3 - ^Je-^-^CNC- (Bl) 

This proves that modes with different momenta are or- 
thogonal. However, for each value of k there are 8 solu- 
tions uj of the dispersion relation (27). Thus, there are 
modes with different frequencies wi ^ UJ2 which share the 



same momentum k„ 



k a2 . We shall show that modes 



with different frequencies uii and 0J2 are also orthogonal. 
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We start from the identity 

C lf v [«0 - ^nk)]v^ = o, 



(B2) 



computed for two modes with the same momentum k ai = 
k a2 = k. Using the relations r/ 2 = —Is, tyrj = /§, and 

(ir,JC V£?n C 2 2 - r,(i V K ) = 0. (B3) 

Since both modes are solutions of Eq. (23), we obtain 

(ui-<»)V£lriV£=0, (B4) 

which shows that V"^ r/V" 2 2 vanishes if uj\ ^ ui2- Since, 
at fixed k ui can take only a finite set of values, it is 
possible to write 



yait„y«2 _<? v ai *nV a2 . 
Putting this expression in Eq. (Bl) 



(B5) 



(C, 1 , Vg ) = 2n6(k a3 - k ai ) 5 U2U1 -V^ ry U« 2 (B6) 



and noting that 

^(^a 2 ^ai ) &L)2U>\ 



dk 



k=k a 



8(u 2 -ui)6 ka2kcii , (B7) 



the combination of Dirac and Kroenecker <5's appearing 
in Eq. (B6) can be rewritten as 



dtu 
~dk 



5(lj 2 - wi) <5 Q2Ql , 



(B8) 



k=k a 



where, as said, a labels the 8 solutions with different 
momentum, sharing the same comoving frequency oj. 

In conclusion, we proved that eigenmodes with differ- 
ent uj and k a are orthogonal: 



did 

dk 



k=k a 



(B9) 

where the last term of this expression is easily computed 
from Eq. (30): 



\Q 



a\2 



2tt 



Eit 



47r/3i7 2 (w + vk) 



(BIO) 

where, for the sake of conciseness, the index a is omitted 
on the right hand side. By using the Lorentz transfor- 
mation (28): 



i v a > n v a = 

UJ I UJ 
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-vK + n 



i + y ^ 

(Bll 



Note that the scalar product is not positive definite. 
Thus, we must characterize under what conditions a 
mode has indeed positive or negative norm. 

From the dispersion relation in the glass rest 
frame (29), the group velocity V g in the laboratory frame 
is given by 



v- 1 = 

y 9 



dK n dK 2 



dn k dn 2 

so that 



47T/3, 



- (i-(i 2 /n 2 ) 2 

(B12) 
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2tt V„ 



vV 



(B13) 



Since there is no absorption, the group velocity \V g \ must 
be smaller than the speed of light c [32]. Furthermore 
v < c, so that the parenthesis in Eq. (B13) is always pos- 
itive. Furthermore, the sign of V g is positive (negative) if 
S7 and K have the same (opposite) sign [see Eq. (B12)]. 
As a consequence, the sign of the above scalar product 
is always equal to the sign of f2, i.e, modes with positive 
(negative) rest frame frequency fl have positive (nega- 
tive) norm. 

Finally, the normalization constant C", appearing in 
Eq. (30) is fixed by imposing 



(B14) 



for both positive and negative norm modes, and by com- 
paring the above condition with Eqs. (B9) and (B13): 
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Using the relativistic composition of velocities 
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1 - VVg/c 2 



Eq. (B13) simplifies to 
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(B16) 



(B17) 



Replacing V g from Eq. (B12) 
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c 2 7 



(B18) 

and going back to the comoving frame frequency u> and 
momentum fc„: 



\C a Y = h 



c k n 
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